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ABSTRACT 

c-j \ The extended TeV gamma-ray source HESS J1023-575 is coincident with the 

^| massive, young stellar cluster Westerlund 2 (Wd2) and its surrounding HII region 

RCW 49. On the basis of an analysis of the CO emission and 21 cm absorption 
along the line of sight to Wd2, it is argued that this cluster, and by assumption 
the TeV source as well, must be associated with a giant molecular cloud in the far 
side of the Carina arm with a mass of 7.5 x 10 5 M Q . Analysis of the spatial and 
velocity structure of the cloud reveals clear evidence of interaction with Wd2. 
The cloud's kinematic distance of 6.0 ± 1.0 kpc is shown to be consistent with 
distances inferred both from the radius-linewidth relation of molecular clouds 
and from the foreground gas column derived from 230 X-ray sources in Wd2. 
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Subject headings: open clusters and associations: individual (Westerlund 2) - 
HII regions: individual (RCW 49) — ISM: individual (HESS J1023-575) — ISM: 
molecules 



1. Introduction 

In just the past few years, a new generation of imaging atmospheric Cherenkov tele- 
scopes, led by the HESS array, have increased the number of confirmed ve ry-high-energy 



(TeV) gamma-ray sources in the Galaxy from just a few (e.g., IWeekesI 120031 ) to more than 
two dozen. Such observations provide insights into where and how highly relativistic parti- 
cles are accelerated and may ultimately identify the main sources of Galactic cosmic rays. 
Some of the Galactic TeV sources can be associated with well-known supernova remnants, 
long considered the most likely cosmic-ray accelerators, and others with pulsars and X-ray 
binaries, while some cannot as yet be associated with any known objects. 
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Unique among the TeV sources so far identified is HESS J1023-575 (lAharonian et al. 



20071 ). an extended source that lies toward Westerlund 2 (hereafter Wd2), a massive stellar 
cluster containing at least a dozen early O stars and two Wolf-Rayet stars. TeV sources and 
such massive clusters are both so rare in the Galaxy and their coincidence in this case so 
close (~ 6') that an association between the two is almost bey ond doubt. The TeV s ource 
J2032+4130 may be associated with a cluster as well (Cyg OB2; Aharonian et al.lbood ). but 
the Cygnus re gion is so complicated that th is association is significantly less secure than in the 
present case ( jButt et al.l 12003c iButtl 120071 ). Various models h ave already been proposed for 



the TeV emission in a ctive stellar regions such as Wd2 (e.g., iBykovl l200ll ; lAharonian et al. 
20071 ; I Bednarek! 120071 ) . most of them involving turbulent particle acceleration driven by 
powerful stellar winds, multiple supernovae, or both. 

A crucial measure for distinguishing among models is the distance to the TeV source, 
which hereafter will be assumed to be the same as that to Wd2. The distance to Wd2, 
however, has a long history of uncertainty, and estimates are not converging with time. Two 



of the most recent distance determinations are among the most discrepant: lAscenso et al. 



(120071 ) used deep near-infrared imaging and ph otometry along wit h pre-main sequence model 
tracks to infer a distance of 2.8 kpc, whereas iRauw et al.l (120071 ) used spectro-photometric 
measurements of 12 O stars in the cluster to estimate a distance of 8.3 ± 1.6 kpc. Most pre- 
vious estimates fall bet ween these values . Thor ough reviews of past d i stance determinations 
were recently given by IChurchwell et al.l (120041 ) and lAharonian et al.l (120071 ) . 



Giant molecular clouds (GMCs) have long provided crucial information about distances 
for a wide range of Galactic objects. GMCs contain a large fraction of the Galactic gas and 
dust, they are generally well defined, are known to be associated with, and often interact 
with, all the products of star formation, and, most importantly, their velocities can readily 
be measured with high precision with CO observations. Here classic techniques in Galactic 
radio astronomy are used in conjunction with CO and 21 cm surveys to argue that Wd2 
and its TeV source are almost certainly associated with a very large GMC lying in the far 
side of the Carina spiral arm at a distance of 6.0 ±1.0 kpc. This view will be supported by 
examination of 21 cm absorption, the foreground extinction, the radius-linewidth relation 
for GMCs, and evidence of interaction between Wd2 and the cloud. 



2. Associated Molecular Cloud 

A reason often cited for difficulty in determining the distance of Wd2 is its direction 
toward the tangent of the Carina spiral arm. However, as Figure 1 shows, Wd2 is actually 
displaced several degrees inward from the tangent, in a direction where the near and far sides 
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of the arm are well separated in velocity from each other and from material near the solar 
circle at km s _1 . Figure 1 also shows that the most conspicuous molecular concentration 
along the line of sight to Wd2 lies squarely on the locus of the far Carina arm at a velocity 
of ~ 11 km s -1 . Here this molecular cloud is called MC8, because it corresponds to cloud 
No. 8 in the Carina arm cloud catalogue of Grabelsky et al. (1988). 

An association between Wd2 and MC8 is strongly implied by an examination of 21 cm 
absorption against the bright radio continuum of RCW 49, the large HI T regio n surrounding 



Wd2. Following the technique described by iMcClure- Griffiths et al.l (120011 ). the average 



21 cm absorption spectrum over the central region of RCW 49 was determined and in Figure 2 
is compared with the CO emission averaged over roughly the same area. It is clear from this 
comparison that all the clouds in the near side of the Carina arm, at negative velocities, as 
well as the cloud at +4 km s _1 , must lie in front of RCW 49, since corresponding 21 cm 
absorption features are detected. Indeed, RCW 49 must lie between the foreground cloud 
F, which produces strong 21 cm absorption, and the background cloud B, which produces 
none. A likely systemic velocity for RCW 49 then is halfway between the velocities of these 
two clouds, which is very close to that of MC8. 

An obvious question raised by Figure 2 is why the strong emission from MC8 itself is 
not evident in the CO spectrum near 11 km s _1 . The answer is provided in Figure 3, which 
displays the emission structure of the cloud both spatially and as a function of latitude and 
velocity. Figure 3b shows that the velocity of MC8 is fairly constant at ~ 11 km s _1 , with 
a gradient of just a few km s _1 from bottom to top. In the general direction of RCW 49, 
however, molecular gas is missing near the systemic velocity of the cloud, and instead clumps 
are seen lying at higher and lower velocities and apparently connected to the main cloud 
by weaker emission. These clumps produce the foreground CO line labeled F in Figure 2 
and the background CO line labeled B. Figure 3b suggests that F and B may arise from 
velocity-perturbed clumps in the same cloud, MC8. 

The energy required to accelerate these clum ps is very modest as compared with that 
available from the known O and W-R stars in Wd2 (IRauw et aU20070 . Assuming N(H 2 ) /W^o 
= 1.8 x 10 20 cm -2 K" 1 km" 1 s ( Dame. Hartmann. Thaddeusl 2001 ) and a distance of 6.0 kpc 



(see below), each clump has a mass of a few 10 4 M , a velocity of ~ 6 km s _1 with respect 
to the main cloud, and a kinetic energy of ~ 2 x 10 49 ergs. This value is ~ 3 orders of 
magnitude less than the wind output of the W-R and O stars over the cluster's estimated 
lifetime of ~ 2 Myr, and one or more past SNe, now extinct, may have contributed energy 
as well. 



The case for identifying clump B as a velocity-perturbed fragment of MC8 is somewhat 
stronger than that for clump F. As Figure 1 shows, emission at the velocity of clump F (~ 4 
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km s _1 ) is not uncommon in the general direction of Wd2, so that a chance alignment of 
F with Wd2 is possible. However, emission peaking above 14 km s _1 is rare in this general 
direction and is seen only toward Wd2 (directly below the Wd2 arrow in Fig. 1). Further, 
the fact that Wd2 and its HII region are optically visible suggests that they lie on the near 
side of MC8, where a general expansion away from Wd2 would lead primarily to a positive 
perturbation of the gas velocities, as is observed for clump B. This molecular clump is the 
one most likely to be closest to Wd2 and the TeV source, p erhaps as close as the core radius 



of the HII region, which is ~ 5' (IWhiteoak fc Uchidal 119971 ). or 9 pc at 6 kpc. Its molecular 



mass is 7 x 10 4 M , mean radius 23 pc, and mean density 30 H2 cm 3 . 

It's worth noting that the proposed association between Wd2 and MC8 is not contingent 
on both F and B being velocity-perturbed fragments of this cloud. The mere appearance of 
a perturbed spatial and velocity structure for MC8 in the direction of Wd2, coupled with the 
tight velocity constraints imposed by the 21 cm absorption spectrum, provide strong evidence 
for association. CO observations with much higher angular resolution will be required to 
determine the true nature of clumps F and B and to investigate in detail the interaction 
between Wd2 and MC8. 

Calculating the H2 mass as above, the total H2 mass of MC8 is found to be 7.5 x 10 



Mm. This is more than twice the mass of the entire Orion molecular complex (IWilson et al. 



2005) and is similar to that of the most massive complexes in the inner Galaxy (IDame et al. 



19861 ). It is not surprising that a cluster as rich in early- type stars as Wd2 is found in one 



of the Galaxy's more massive molecular complexes. 



Distance 



Determining the distance to MC8 is in principle straightforward, because the positive 
velocity of MC8 places it beyond the solar circle, where, on the assumption of pure circular 
rotation, there is a monotonic relationship between velocity and distance. However, the 
Galactic rotation curve is more uncertain beyond the solar circle than within, and one 
must consider, in addition, the uncertainty imposed by random and systematic non-circular 
motions. Given that the well-determined rotation curve inside the solar circle is flat to within 



a few km s _1 kpc _1 (jBrand fc Blitall993l ) and that the rotatio n curves of external spirals are 
generally flat or slowly rising to well beyond the solar radius (IBinney fc Merrifieldlll998l ). it 
is very unlikely that the rotation curve between the solar circle and MC8 has a significant 
slope. Us ing more than 200 H II regions and reflection nebulae with spectro-photometric 
distances, iBrand fc Blitzl (119931 ) derived a slope of ~ 1 km s _1 kpc~ 1 for the rotation curve 
beyond the solar circle. Use of this curve with R = 8.5 kpc and a velocity of 11 km s _1 for 
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MC8 implies a distance of 6.0 kpc. Changing the slope of the rotation curve of the outer 
Galaxy by an improbably large amount with respect to the Brand & Blitz curve, ±5 km 
s _1 kpc _1 , changes the derived distance by only ±0.35 kpc. 

In considering the effects of random and systematic noncircular motions, it is important 
to keep in mind that the Galaxy appears overall to be a very orderly system. The smooth 
and well-defined loci of terminal velocities with longitude in the first and fourth quadrants , 
and their agreement between quadrants to better than 10 km s~ 1 (jBlitz fc Spergell Il99ll ). 
tightly constrain the amplitudes of noncircular motions. Random jitter about these smooth 
curves implies a one- dimensional velocity dispersion of ~ 4.2 km s _1 for molecular clouds 
(jStark fc Brand! 119891 ). and broad "bumps" on the curves at the spiral arm tangents imply 
an am plitude of ~ 7 km s _1 for density-wave streaming along the arms (jBinney &: Merrifield 
19981 ). Larger velocity perturbations may occur perpendicular to the arms because of spiral 
shocks (jRobertslll969l ). but at MC8 a shock perturbation would be directed at a large angle 
with respect to our line of sight, ~ 74° assuming an arm inclina tion of 1 ° (IGrabelsky et al. 



19881 ). The present case is in contrast to the seminal study by iRobertd (119721 ) of shocks in 
the Perseus arm in the second quadrant, where the proposed shocks of amplitude ~ 20 km 
s _1 are directed almost exactly toward the Sun. A similar shock at the position of MC8 in 
the Carina arm would perturb its velocity along the line of sight by at most 5 km s _1 . The 
actual velocity perturbations due to streaming and shocks depend on the exact location of 
MC8 with respect to the potential minimum of the Carina arm, which is unknown, but the 
perturbations should be smaller than the full streaming and shock amplitudes. Therefore, 
random motions, streaming, and shocks together imply a velocity uncertainty for the cloud 
of ~ 8 km s~ 1 and a corresponding distance uncertainty of ±0.8 kpc; this increases to ±0.9 
kpc when the uncertainty of the rotation curve in the outer Galaxy is included. Thus one 
can conclude that MC8 along with Wd2, RCW 49, and HESS J1023-575 all lie at a distance 
of 6.0 kpc with a conservative uncertainty of ±1.0 kpc. 

A rough independent check on the distance to MC8 can be ob tained from the radius- 
linewidth relation derived for Galactic GMCs by iDame et al.l (119861 ). Following this study, 
the cloud's angular radius is calculated as (A/tt) 1 ^ 2 , where A is the projected area of the 
cloud as in Figure 3a. Excluding the perturbed section of the cloud toward Wd2 (see Fig. 
3b), one obtains a composite cloud linewidth of 8.3 km s _1 and, at 6.0 kpc, a linear radius of 
50 pc. As Figure 4 shows, a distance of 6.0 kpc places Wd2, perhaps somewhat fortuitously, 
almost exactly on the fit line. Distances in the range 4-8 kpc are allowed by the fairly large 
scatter. 



An independent check on the distance to Wd2 can be obtained by comparing its fore- 
ground gas column as derived from X-ray absorption to that derived from the gas tracers. 
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Very recently, iTsujimoto et al.l (120071 ) identified 468 X-ray sources in Wd2 with Chandra and 
fit each of the stronger 230 of these with an optically thin thermal plasma (APEC) model 
to derive foreground N(H). The sources are spread over a cluster core radius of 6' to 7', a 
good match to the 8.4' resolution of the CO data used here. The distribution of the derived 
N(H) is strongly peaked at 1.5 ± 0.5 x 10 22 cm -2 (see their Fig. 5). As discussed above, the 
21 cm absorption spectrum (Fig. 3) demonstrates that all the CO and 21 cm emission at 
v < 10 km s _1 must lie in front of Wd2. Integrating the nearest CO spectrum to Wd2 over 
velocities of less than 10 km s _1 yields N(H 2 ) = 0.44 x 10 22 cm " 2 . Sim ilarly, integrating the 
nearest 21 cm spectrum from the HI survey of iKalberla et al.l (120051 ) . one obtains N(H) = 
0.56 x 10 22 cm -2 . A mean molecular weight per H of 1.36 implies a total foreground N(H) 
of 1.9 x 10 22 cm -2 , in reasonable agreement with the value derived from the X-rays. 

Although it is not possible here to discuss all previous dista nce estimates for Wd 2, it is 
worth noting that the lowest estimate frequently cited (2.2 kpc; iBrand fc Blita 119931) actu - 
ally derives from an unpublished spectro-photometric study in the thesis by iBrandl (119861 ). 
Another group of low estimates, of order 4 kpc, are kinematic distances based on radio 
recombination- line measurements of the velocity of RCW 49, which are ~ km s _1 (e.g., 



Wilson et al.lll970l ; iBenaglia et al.l 120051 ) . It has long been known that optically visible HII 
regions such as RCW 49 tend to be blueshi f ted w ith respect to their parent clouds in a manner 
consistent with the blister model of Israel Jl978h . Given the evidence of interaction between 
RCW 49 and MC8 presented above, it seems likely that the relatively low-mass ionized gas 
will be perturbed toward us by an amount larger than the 6 km s _1 perturbation of molecular 
clump B away from us. Such a blister flow would shift lower both the recombination-line 
velocity of RCW 49 and its derived kinematic distance. 



4. Summary 

The HESS collaboration has recently identified a new type of very-high-energy gamma- 
ray source, the massive young stellar cluster Wd2. Attempts to understand the extended 
TeV emission from this complicated region have so far been hampered by a long history of 
discrepant distance estimates for Wd2, with recent values ranging from 2 to 8 kpc. 

Comparison of a CO emission spectrum toward Wd2 with a 21 cm absorption spectrum 
measured against its surrounding HII region RCW 49 indicates that Wd2 must lie behind a 
molecular concentration at +4 km s _1 and in front of one at +16 km s -1 . The only substantial 
molecular concentration in this velocity range is also the most massive cloud along the entire 
line of sight, a cloud here called MC8. The Wd2 cluster is so massive and young that it would 
be expected to lie very near a large reservoir of molecular gas such as MC8. In the direction 
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of Wd2, MC8 appears to be substantially disrupted both spatially and in velocity One or 
both of the molecular concentrations at +4 and +16 km s _1 might be velocity-perturbed 
fragments of MC8. 

At a velocity of 11 km s _1 , MC8 lies along the longitude- velocity locus of the far Carina 
arm, at a kinematic distance of 6.0 kpc. Uncertainty about the rotation curve in the outer 
Galaxy and possible deviations from pure circular motion due to density-wave streaming, 
spiral shocks, and random cloud motions cumulatively imposed an uncertainty of ±1.0 kpc 
on the kinematic distance. At 6 kpc, the radius and linewidth of MC8 are in good agreement 
with the radius-linewidth relation derived for large inner-Galaxy clouds, and the foreground 
extinction derived from the gas tracers is in good agreement with that derived from 230 
X-ray sources in Wd2. 

I am grateful to Naomi McClure- Griffiths for data and advice on the Southern Galactic 
Plane HI survey, Paula Benaglia for data and useful discussions on radio recombination lines, 
and Yousaf Butt for bringing the discovery of HESS J1023-575 to my attention. 
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Fig. 1. — Longitude- velocity map of CO emission integ rated over a strip ~ 4° wide in l atitud e 
centered on the Galactic plane. The data are from iDame. Hartmann. Thaddeusl (120011 ). 
smoothed slightly to a resolution of 2 km s _1 in velocity and 12' in longitude to emphasize 
the large-scale structure of the Carina arm, which appears as a large loop running from 
upper left to lower left. Intensities are in units of log(K-arcdeg), ranging from 1 (dark blue) 
to 5 (white). 
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Fig. 2. — Comparison of CO emission ( IGrabelsky et al.1 Il988l ) and HI absorption 
(IMcClure- Griffiths et al.ll2005l ) toward the giant HII region RCW 49, which surrounds Wd2. 
The 21 cm spectra from the brightest part of the continuum emission (I = 284.3° to 284.35°, 
b = —0.38° to 0.32°) were averaged as the "on" and spectra offset 5' to 10' from the con- 
tinuum peak were averaged as the "off". The optical depth spectrum was obtained by 
subtracting the "off" from the "on" and dividing by the off-line continuum level of the "on" . 
The CO spectra are averaged within ~ 10' of the continuum peak. The letter F marks the 
highest- velocity line from molecular gas in front of RCW 49, and B marks the lowest- velocity 
line from molecular gas behind RCW 49. 
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Fig. 3. — (a) Color: CO intensity integrated from to 20 km s _1 in LSR velocit y. Contours: 
Cont inuum emission at 843 MHz from the Molonglo Galactic Plane Survey (IGreen et al. 
19991 ): contours are in units of log(mJy beam -1 ), from 1.5 to 3.6 in steps of 0.3. The 
red dotted circle marks the position and FWHM size of HESS J1023-575. The shell-like 
continuum source centered near 283.9-0.9 is RCW 48, which may be related to RCW 49 and 
MC8. The CO source at 285.25+0 is probably unrelated to MC8. (b) Latitude- velocity map 
of CO intensity integrated over the longitude range of MC8, I = 284° to 285°. The letters F 
and B mark possible fragments of MC8 expanding away from RCW 49; as Fig. 2 shows, F 
is seen in absorption against RCW 49, but B is not. 
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Fig. 4. — The radius-linewidth relation for giant m olecular clouds in the first Galactic quad- 
rant (adapted from Fig. 6 in iDame et al.l (119861 )). The horizontal line marks the locus of 
possible positions of MC8 for assumed distances ranging from 2 to 8 kpc. MC8 fits best 
the radius-linewidth relation at a distance of 6 kpc, the same distance that was derived 
independently from the cloud's velocity. 



